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18.  SU8JECT  TERMS  (Continued) 

HF-scatter  Model 
Scatter-propagation  Model 


SUMMARY 


M  V.  V. 


TKIV ”•  VA1 1 
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Methods  have  been  described  for  evaluating  the  major  effects  on  HF 
sky-wave  signals  of  forward  scatter  from  ionospheric  irregularities. 

The  approach  used  to  obtain  these  results  combines  the  scatter  theory 
that  Flattg,  et  al.1,  developed  for  refracted  ray  paths  and  that  developed 
under  DNA  sponsorship  ’  for  transionospheric  radio-wave  propagation. 

The  effects  addressed  by  the  theory  include  angular  spread/spatial 
decorrelation,  frequency  (Doppler)  spread/time  decorrelation,  and  time 
spread/frequency  decorrelation.  The  latter  two  of  these  effects, 
frequency  spread  and  time  spread,  are  those  conventionally  represented 
by  the  channel-scattering  function.  Angular  spread  is  of  concern  to  the 
HF  system  analyst  for  systems  with  narrow-beam  antennas  or  spatial- 
diversity  capabilities. 

The  results  presented  here  incorporate  a  number  of  simplifications 
thought  to  be  appropriate  in  applying  the  theory  in  a  system  context. 
Improving  approximations  in  this  spirit  (e.g.,  to  provide  a  better 
representation  of  the  vertical  profile  of  ionospheric  electron-density) 
would  in  many  cases  be  straightforward,  but,  at  the  expense  of  addi¬ 
tional  complexity  in  the  algebraic  detail. 

A  number  of  areas  merit  further  attention.  Evaluating  the 
frequency  coherence  bandwidth  (Section  7)  may  well  require  elaboration. 
However,  evaluating  the  theory  in  its  present  form  would  seem  appro¬ 
priate  first.  The  HF  propagation  experiment  currently  being  undertaken 
by  SRI  International  under  DNA  sponsorship  addresses  this  need. 

Comparing  the  effects  observed  experimentally  with  those  predicted  by 
the  present  theory  will  provide  the  best  guide  to  those  areas  in  which 
further  development  will  be  worthwhile. 

A  question  of  recurring  interest  concerns  propagation  by  means  of 
paths  reflected  from  nuclear  plumes.  These  regions  of  enhanced. 


i  i  i 


structured  ionization  are  produced  by  nuclear  detonations  in  the  upper 
atmosphere.  The  theory  presented  here  appears  capable  (with  a  suitable 
redefinition  of  what  is  meant  by  "background  ionosphere")  of  providing 
means  to  assess  some  of  the  characteristics  of  such  modes. 
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SECTION  1 


INTRODUCTION 

1  * 

Over  the  past  decade,  it  has  been  recognized  that  the  scattering 

of  signals  by  a  randomly  irregular  medium  can  be  characterized  by 

integrals  evaluated  along  the  unscattered  ray  paths.  This  work 

can  accommodate  inhomogeneous  random  media,  in  which  the  ray  path  is 

refracted  even  in  the  absence  of  irregularities  and  anisotropic 

irregularity  structures.  The  application  of  these  results,  together 

with  the  model  of  ionospheric  irregularities  developed  by  the  Defense 

2  3 

Nuclear  Agency  (DNA)  for  application  to  satellites,  *  to  the 
propagation  of  ionospherically  refracted  high-frequency  (HF)  radio  waves 
is  described  in  this  final  report  for  Contract  DNA  001  -83-C-0413. 

In  addition  to  developing  procedures  for  calculating  scatter 
effects  reported  here,  this  contract  also  included  tasking  to  support 
the  planning  for  the  HF  Channel  Probe  Experiment  that  was  fielded 
successfully  in  Greenland  by  SRI  International  during  the  Fall  of  1984. 
Developing  the  theory  in  conjunction  with  the  experimental  program  has 
significantly  enhanced  the  utility  of  both  endeavors.  Evaluating  theory 
presented  here  through  the  comparison  with  the  experimental  results 
should  be  substantially  easier  because  of  the  interaction  between  the 
experiment  and  the  modeling  development.  Indeed,  as  a  result  of  this 
interaction,  selecting  the  quantities  to  be  measured  experimentally  and 
specifying  the  particulars  of  their  measurement  were  made  with  close 
attention  to  achieving  meaningful  and  straightforward  comparisons  with 
theory. 


* 

References  are  listed  at  the  end  of  this  report. 


Three  quantities  generally  suffice  to  characterize  the  effects  of 
irregularities  adequately  for  most  system  applications.  These  are  the 
correlation  between  samples  of  the  received  signal  separated  (1)  tempor¬ 
ally,  (2)  spectrally  (i.e.,  in  frequency),  and  (3)  spatially.  Attention 
is  first  given  to  calculating  the  phase-structure  function.  This 
function  is  simply  related  to  the  spatial  decorrelation  of  the  signal 
phase  on  two  paths  separated  by  a  small  displacement  at  the  receiver. 
This  signal  decorrelation  between  two  separated  receiving  antennas 
determines  the  effectiveness  of  spatial  diversity  in  improving  system 
performance.  The  spatial  correlation  function  is  also  directly  related 
to  the  angular  spread  of  the  received  signal.  Although  HF  antennas  are 
not  easily  made  highly  directive,  some  have  been  so  designed,  and  the 
degree  of  angular  spread  relative  to  the  antenna  directivity  remains  an 
important  parameter  in  assessing  HF  system  performance. 

The  methods  described  for  calculating  the  phase-structure  function 
also  pertain  to  calculating  other  scatter  parameters.  As  will  be  shown, 
only  minor  modification  of  the  procedure  is  necessary  to  obtain  a 
description  of  the  temporal  decorrelation  of  the  signal,  given 
reasonable  (with  our  present  knowledge)  simplifying  assumptions. 

Indeed,  the  concept  of  the  phase-structure  function  can  be  generalized 
to  include  separations  in  time  and  frequency  as  well  as  space. 

Calculating  the  phase-structure  function  requires  determining  the 
path  of  the  unscattered  ray.  Determining  the  sensitivity  of  the 
integrated  phase  along  the  path  to  a  displacement  of  the  path  from  its 
nominal  (unscattered)  position  is  also  necessary.  To  facilitate  these 
calculations,  we  have  simplified  the  description  of  the  background 
(i.e.,  unstructured)  ionosphere  somewhat.  This  simplification, 
described  in  Section  2,  is  a  conventional  one  often  employed  in 
ionospheric  radio-wave  studies:  The  ionospheric  electron-density 
distribution  is  assumed  laterally  uniform  (i.e.,  the  ionosphere  is 
assumed  stratified),  and  the  form  of  the  variation  in  electron  densitv 
with  height  is  chosen  such  that  the  refracted  ray  paths  can  be 
analytically  described  (Section  3).  Thus,  the  computational  burden 
imposed  by  incremental  rav-trace  solutions  for  the  ray  paths  between 


specified  points  have  been  avoided.  We  anticipated  as  well  that  our 
needs  would  have  imposed  precision  requirements  on  an  incremental  rav- 
trace  approach  that  substantially  exceeded  those  normally  encountered. 
With  the  analytic  result,  it  is  a  simple  matter,  for  small  displace¬ 
ments,  to  determine  the  displacement  to  the  ray  path  for  any  point  along 
the  path  given  a  specified  displacement  of  an  end  point,  such  as  is 
needed  for  calculating  the  phase-structure  function  (Section  4),  without 
repeating  the  entire  path  calculation. 

Also  described  briefly  in  Section  2  is  the  model  used  to  represent 
the  irregularity  structure  of  the  ionosphere  responsible  for  the  effects 
of  scatter.  This  model,  which  has  been  adapted  from  that  developed  as 
part  of  DNA's  efforts  to  characterize  the  effects  of  scatter  at  higher 
frequencies,  seems  a  reasonable  starting  point  for  describing  the 
ionospheric  structure  important  to  HF.  Elaboration  or  modification  of 
the  irregularity  model  may  prove  desirable  as  experience  with  the  HF 
scatter  model  is  gained.  Because  data  are  lacking  with  which  to  address 
this  question  at  the  present  time,  we  will  not  address  it  here. 


Characterizing  the  signal  decorrelation  in  space,  time,  and 
frequency  is  discussed  in  Sections  5  through  7  after  the  development  in 
Section  4  of  methods  for  calculating  the  phase-structure  function.  We 
find  that,  as  in  the  simpler  case  of  transionospheric  propagation  at 
higher  frequencies,  a  nominal  value  for  the  power-law  index  of  the 
irregularity  spectrum  leads  to  a  simple  form  for  the  spatial  and 
temporal  correlation  functions.  The  decorrelation  of  the  signal  in 
space  and  time  consequently  can  be  simply  characterized  by  decorrelation 
length  (which,  however,  depends  on  direction)  and  decorrelation  time, 
respectively.  The  relationship  of  the  spatial  decorrelation  length  to 
the  angular  extent  of  the  signal,  which  is  a  more  directly  useful 
quantity  in  many  cases,  also  is  described. 


V 


SECTION  2 


IONOSPHERIC  MODEL 

We  have  assumed  that  the  unperturbed  ionospheric  electron  density 
varies  only  with  height.  The  form  of  this  variation  is  considered 
quasi-parabolic,  a  choice  that  is  analytically  convenient  without 
sacrificing  practical  utility.  Specifically,  we  assume  a  spherically 
stratified  profile  of  the  form 

(  N  [l-{  (r-r  )/y  }2(r,  /rf]  ;  r,  <  r  <  r.r  /(r.-y  ) 

1  lv  K  b  ;  J  b—  —  b  m  v  b  m ' 

N  (r)  =  v  (I) 

e 

'  0  ;  elsewhere 

in  an  earth-centered  spherical  polar-coordinate  system  with  r  the  radial 

distance  from  the  earth's  center.  The  layer  peak  is  located  at  r  ,  with 
_  m 

a  peak  density  N  ,  and  its  base  at  r, ;  the  layer's  semithickness  is 
m  d 

v  •  r  -  r,  .  For  this  profile,  ray  paths  can  be  described  analyt- 
•  m  m  b  J  r 

4 

ically  by  a  set  of  equations  giving  the  angular  distance  along  the 
earth's  surface,  0(r),  as  a  piecewise  continuous  function  of  r.  As  was 
noted  in  the  introduction,  the  availability  of  these  equations 
facilitates  calculating  the  scatter  parameters  for  the  HF  channel 
model. 

The  ionospheric  irregularities  are  assumed  anisotropic,  as 

described  by  a  three-dimensional  power-law  spatial  spectrum  of  electron- 

3 

density  fluctuations  of  the  form 


a 


4>(q)  = 


v+1/2 


where  q  is  a  quadratic  form  in  the  wave-number  components  K  ,  K  ,  and 

K  ;  v  is  the  power-law  index;  and  a'  and  b'  are  elongation  factors  that 
z 

describe  the  correlation  lengths  in  the  directions  of  the  second  and 

third  principle  axes  relative  to  the  shortest  in  the  direction  of  the 

first.  Here,  a'  describes  the  elongation  parallel  to  the  earth's 

magnetic  field,  and  b',  that  perpendicular  to  the  geomagnetic  meridian 

plane;  both  are  normalized  to  the  correlation  length  perpendicular  to 

* 

the  field  direction  in  the  meridian  plane. 

The  parameter  Cg  is  a  measure  of  irregularity  strength,  and  q^ 

specifies  the  outer  scale.  The  electron-density  variance,  ,  can  be 

obtained  by  integration  of  $(q)  over  all  q.  The  rms  electron-density 

variation,  o  ,  is  considered  a  fixed  fraction  of  the  mean  electron 
N 

density,  <N  >;  such  a  relationship,  although  not  essential  to  the  model, 
e 

is  a  reasonable  approximation  for  the  propagation  environments  of 

interest.  With  q  and  N  specified,  C  can  be  calculated  from  the 
o  s 

formula 


„  3/2  2  2 v-2  .  ^ 

Hit  o..q  T(v  +  1/2) 
_ N  o _ 

T(v  -  1) 


The  coefficients  of  K  ,  K  and  K  depend  on  the  coordinate  system,  as 
will  be  described  shortly^ 


Figure  1  relates  the  numerical  magnitudes  of  0^  and  for  various 

values  of  the  outer  scale  L  =  L/q  . 

o 


Cs  —  mks  units 


Figure  1.  Conversion  of  rms  measures  to  turbulent  strengths. 

The  relationship  of  the  parameters  defined  here  to  those  employed 
by  Wittwer^  to  model  structured  ionization  is  summarized  in  Table  I. 
The  two  models  are  formally  identical. 


Table  1.  Notation  for  propagation  modeling 


Natural 

2 

Ionosphere 

Wittwer 

% 

L-1 

s 

a' 

L  /L 
t  s 

b' 

L  /L 
r  s 

<5 

6 

G 

q  L 

O  Z 

* 

^BP 

<t> 

aqo/G 

,  ,2.2  2  .  2  . 

L  =  L  (a  cos  *  +  sin 
u  s 

0 

9 

z  sec 

K 

9 

<ZT  '  2l,Zl 

ZC 

f  zT 

V 

n  =  v  +  1/2 

AN 

e 

-(2 v+1 ) 
q 

.  -  2n 

15  AN  q 

e 

Comment 


Outer-scale  wave  number 

Axial  ratio  along  field 

Transverse  axial  ratio 

Orientation  of  transverse 
irregularity  axis 

Geometric  enhancement 
factor 

Briggs-Parkin  angle 
Br iggs-Parkin  factor 
Zenith  angle 
Magnetic  azimuth  angle 

"Reduced"  propagation 
distance 

Spectral  index  parameter 

Form  of  three-dimensional 
spectral  density  function 


The  Br iggs-Parkin  angle — the  angle  between  the  line  of  sight  and  the 
magnetic  field — is  given  in  terms  of  if,  9,  and  <f  as 


ros  if  =  cos  <1  sin  0  cos  <i>  +  sin  T  cos  A 


SECTION  3 

UNPERTURBED  RAY  PATP 


aiLUi 


>7>.l 


In  order  to  compute  the  effects  of  structure,  a  description  of  the 
unperturbed  ray  path  is  needed.  As  will  be  described  shortly,  these 
computations  involve,  among  other  things,  determining  the  extent  to 
which  the  path  is  shifted  by  a  displacement  of  its  end  points.  Calcu¬ 
lating  this  displacement  is  made  easier  by  the  form  selected  for  the 
unperturbed  electron-density  profile,  Eq.  (1);  this  form  allows  the 
unperturbed  path  to  be  specified  as  an  explicit  function  relating  the 
signal  frequency,  profile  parameters,  and  path  coordinates.  In  this 

development,  we  have  attempted  to  retain  the  notation  of  Croft  and 
4 

Hoogasian  for  the  profile  and  path  parameters. 

Snell's  law  for  a  spherically  stratified  medium  (also  known  as 
Bourger's  rule)  specifies  generally  that  the  local  elevation  angle  of 
the  ra',  3,  satisfies 


B(r,  3  )  =  arc  cos  [r  cos  3  /(rjj(r))]  , 


(A) 


where  3  is  the  ray  elevation  angle  at  the  transmitter.  The 
o 

(unperturbed)  mean  refractive  index,  y,  is  given  by 


-2 ,  .  .  e2Vr) 

V  (r)  =  1 - 2 - 

e  u)  m 
o  e 


(5) 


For  the  region  between  the  transmitter  at  r  =  (r  ,  0,  0)  and  the 

t  o 

base  of  the  ionosphere  at  r^,  this  relationship  yields 


9(r,  3  )  =  3(r,  3  )  -  3 
o  o  o 


r  <  r  <  r, 
o—  —  b 


(6) 


r.v.M 

mm 

-•  *.* 


*1 
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while  within  the  ionosphere, 


r  cos  6  , 

9(r,  0  )  =  0,(0  )  -  0 - - —  Un[n(r)  sin  0(r,  0  ) 

o  b  o  o  cl/2(0  )  l 

o 

+  c1/2(0  )/r  +  b/(2c1/2(0  ))] 
o  v  o 

-  £n[sin  0b(0o)  +  cl/2(0o)/rb  +  b/ f 2 c 1 /2 ( ) ] }  .  rfe  <  r  <  rt 

with 


and 


0  (0  )  =  0(r.  ,  0  ) 
bo  bo 


(8) 


c(  0  )  =  (r,r  /(Fy  ))2-  r2  cos20 
o  v  b  m  •'ra  ‘  o  o 


(9) 


where 


F 


(10) 


is  the  ratio  of  the  wave  angular  frequency,  w,  to  the  critical  frequency 


2 

U) 

c 


'  e2N 


m 


e  m 
o  e 


(11) 


The  parameter  b,  not  to  be  confused  with  the  irregularity  anisotropy 
parameter  b'.  Is  given  by 


4 


b  =  -  2—  ( r  r  /(Fy  ) ) 
r  v  b  m  ■'m 


On  the  descending  leg  of  the  ray  path,  the  corresponding 
expressions  for  9  are 


r  cos  8 


9(r,  6  ) 

o 


-  MB  )  -  0  - 


-V  \  M  I/O 

bo  o  c1/2(8  ) 

o 


£n  [wi "  a] 

o 


c1/2<8> 
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(15) 


again  is  not  to  be  confused  with  the  irregularity  anisotropy  factor  a 
The  height,  r  ,  of  the  rav  apex  is  given  bv 


SECTION  4 


PHASE-STRUCTURE  FUNCTION 


The  mutual  coherence  function  of  a  randomly  scattered  wavefield 
describes  the  wavefield  decorrelation  caused  by  the  scatter.  Spatial, 
temporal,  and  spectral  decorrelation  of  the  signal  can  all  be  important 
to  a  communication  system  and,  therefore,  need  to  be  described  in  the 
channel  model.  Spatial  decorrelation  is  conveniently  expressed  in  terms 
of  the  phase-structure  function,  D(p,p"),  which  is  defined  as  the  mean 
square  difference  in  the  signal  phase  perturbation,  6<J>,  observed  at  two 
points,  p  and  p' : 


D(  A)  -  <[  6<t>(P )  -  6<Kp')]2>  ,  (17) 

where  A  =  |p  -  p'\  is  the  spatial  separation  between  the  points. 

Temporal  coherence  of  the  wavefield  is  also  usefully  discussed  in  terms 
of  a  similar  quantity  if  the  temporal  variation  of  the  propagation 
medium  is  simply  equivalent  to  a  uniform  drift  of  the  irregularities. 

We,  therefore,  first  consider  the  evaluation  of  the  phase-structure 
function.  In  following  sections,  we  address  using  it  and  related 
quantities  to  characterize  the  signal  decorrelation. 

2  5 

The  phase-structure  function  can  be  approximately  evaluated  ’  as  a 
sum  over  contributions  from  points  along  the  unscattered  ray  path 

D(AR)  =  /  ds  d [ r ;  A(r)]  (18) 

path 

where  A  =  A(r  )  =  p  -  p'  is  some  spatial  displacement  at  the  receiver, 

K 

located  at  r  ,  and  A(r)  is  the  corresponding  displacement  at  some  point, 

”  R 

r  =  r(s),  on  the  ray  path.  The  relationship  of  A(r)  to  A  is  discussed 


further  below. 


Because  0  is  given  explicitly  as  a  function  of  r  by  the  ray-path 
equations,  it  is  straightforward  to  express  the  phase-structure 
function,  D,  as  an  integral  over  r.  The  increment  of  path  length  is 
given  by 


sin  0(r) 


and  D(A  )  can  be  written 


°<*R>  ■ 

rb 

r rb ,  d{ r;  A(r)] 

+  J  dr  Tin  0(Y) 


if  we  assume  that  the  path  segments  between  the  earth's  surface  and  the 
base  of  the  ionosphere  are  free  of  irregularities. 

The  incremental  contribution  to  the  phase-structure  function  from 
each  differential  segment  of  the  path  can  be  viewed  as  equivalent  to 
that  from  a  phase  screen  representing  the  effects  of  the  irregularities 
located  within  a  layer  of  thickness  dL  =  dr/tanB.  The  general  form  of 
the  phase-structure  function  for  such  a  layer  has  been  calculated  by 
Rino^;  formally,  we  equate  this  result  to  dir;  A(r)]ds  to  obtain 


d[r;  A(r) ]  = 


r2X2G(r)C  (r)T(v  -  1/2) 
e  s 


2*r(v  +  1/2) 


-•(¥) 
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K 

Vl/2  ° 
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2  v  -  1 


where  X  is  the  signal  wavelength  and  r  is  the  classical  electron 

e 

radius.  t*ie  m°dified  Bessel  function  of  imaginary  arguments, 


-  C A,  -BA .A  +  AA 
2  _  4>  <frr  r 

"  2 

AC  -  B  /4 


A  =  C. 


(24a) 


B  =  2  (C.,  -  tan|e|C. ,) 


(24b) 


C  =  C22  -  2C23  tan  |B|  +  C33  tan  B  (24c) 


providing  the  dependence  of  the  scatter  effects  on  the  path  orientation 
relative  to  the  geomagnetic  field.  This  dependence  is  introduced  by  the 
anisotropy  of  the  irregularities. 


Equation  (21)  can  be  approximated  well  by 


dir;  A(r)  ]  =  G(r)C^.  |y  |inln(2v  l’2)  , 


where 


.  2  2 

X  reCs  ,  2  9  .  1 1 > 

2*  (V  '  2  V  +  t1  ' 


An  equivalent  approximation  has  been  derived  by  Wittwer.  ’ 


The  C. .  are  the  elements  of  a  transformation  matrix  for 

ij 

normalization  of  a  quadratic  form  that  appears  in  the  anisotropy  model. 
For  our  geometry. 


(27a) 


(sin2*?  +  a'2  cos24')  sin26  +  b'2  cos2S 


2  2  2 
cos  f  +  a'  sin  f 


(27b) 


(sin2*}1  +  a"2  cos2^)  cos25  +  b'2  sin25 


(27c) 


=  (1  -  a'  )  sin  V  cos  41  sin  6 


(27d) 


=  (sin2^  +  a"2  cos2  31  -  b'2)  sin  6  cos  6 


(27e) 


=  (1  -  a'  )  sin  41  cos  V  cos  6  , 


(  2  7  f  ) 


where  V  is  the  dip  angle  of  the  geomagnetic  field,  and  6  is  the  angle 
between  the  geomagnetic  meridian  plane  and  that  containing  the  path.  A 
positive  6  is  taken  to  rotate  0  into  <j>  as  the  path  plane  is  rotated 
about  i^  into  the  geomagnetic  meridian  plane. 

Given  a  displacement  A  =  [A^,  Ar>  Ag)  at  the  receiver,  the  A^  and 

A  are  functions  of  the  location  on  the  ray  path.  First,  the 
T  R  |  11 

displacement  A  is  reduced  to  that  in  the  (r,  4>)  plane,  A  =  (a.,  A  1 

<P  r 

through  the  relationship 


A1  =  (A*  AR  +  tan  0  AR) 
v  <j)  r  o  0' 


The  displacement  A^(r)  in  the  horizontal  transverse  coordinate  is  simply 
proportional  to  0: 


.  ,  .  0(r)  .1  _  0(r )  ,R 

*  ‘  V  *  'V  *  ' 


The  displacement  Ar(r)  in  the  vertical  coordinate  is  less  direct. 

Given  0(r,  8  ),  we  have 
o 

d8  =  dr  +  ||  d8 
9r  30  o 
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(30) 


for  the  displacement  in  0,  given  incremental  changes  in  r  and  the  take 

off  angle  0  .  For  a  fixed  0,  d0  =  0,  and  r  is  related  to  0  by 
o  o 

,  (30/30  ) 

i£_  =  -  _ _ (31) 
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For  small  A^,  then,  the  change  in  0q  induced  by  the  displacement  at 

0  causes  a  displacement  A  (r)  at  9(r)  given  by 
R  r 
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Calculation  of  the  partial  derivatives  of  9  from  the  expressions  given 
above  for  9(r,  0  )  is  messy,  but  straightforward. 

Only  the  practical  difficulties  involved  in  evaluation  of  the 
integrals  for  D(A  )  remain  to  be  considered.  At  the  path  apex,  6=0. 
Consequently,  we  have  to  deal  with  an  integrable  singularity  resulting 
from  the  sin  0  in  the  denominator.  This  is  a  purely  geometric  factor 
resulting  from  the  projection  of  the  path  increment,  ds,  onto  the 
vertical-axis  integration  increment,  dr.  The  offending  factor  can  be 
cast  in  the  form 


1  u(r)r _ 

sin0( r )  ",  2  A  .1/2 

(ar  +  br  +  c) 


showing  that,  given  the  absence  of  any  other  singularities  at  r 
introduced  by  the  remainder  of  the  integrand,  the  singularity  can  be 


integrated  (cf . ,  the  integral  forms  for  other  path  characteristics 

4 

addressed  by  Croft  and  Hoogasian  ).  Although  algebraically  messy, 
expansion  of  the  remainder  of  the  integrand  in  a  power  series  about  r 
provides  a  straightforward  approach  to  dealing  with  the  singularity. 
This  expansion  is  employed  over  only  a  limited  portion  of  the 
integration  range,  with  straightforward  numerical  integration  employed 
over  the  remainder,  to  gain  efficiency. 


SECTION  5 


SPATIAL  DECORRELATION 


Spatial  decorrelation  of  the  received  signal  is  described  directly 
by  the  phase  structure  function,  D(A),  defined  in  Eq.  (17).  Namely,  the 
spatial  autocorrelation  function  for  the  complex  signal  envelope  is 
given  by 


-  ^D(A) 

R( A)  =  e  .  (35) 

Although  the  spatial  correlation  function  is  of  immediate  interest  in 
the  assessment  of  spatial-diversity  systems,  the  angular  power  spectrum 
of  the  received  signal  is  usually  more  directly  useful.  The  two 
quantities  constitute  a  Fourier-transform  pair  when  the  angular  spectrum 
is  expressed  in  wave-number  (k)  space.  Thus,  one  can  readily  be  derived 
from  the  other,  and  both  contain  the  same  information. 

For  sufficiently  small  spatial  displacements,  the  relationship 
between  the  displacement  of  the  path  end  point  and  that  of  other  points 
along  the  path  is  linear,  as  given  by  Eqs.  (29)  and  (33).  Equation  (29) 
remains  linear  regardless.  In  the  linear  case,  the  dependence  of  D(A) 
on  A  is  the  same  as  that  given  by  Eq.  (25)  for  d[r;A(r)],  the 
incremental  contributions  to  D  from  different  portions  of  the  path. 

Thus,  given  the  commonly  assumed  value  of  1.5  for  v,  the  power-law  index 
of  the  ionospheric  irregularities,  D  is  a  quadratic  function  of  A,  and 
the  autocorrelation  function  is,  from  Eq.  (35),  Gaussian  in  form. 

The  Fourier  transforms  relating  angular  power  spectrum  and  spatial 
correlation  function  can  readily  be  evaluated  under  these  circumstances, 
yielding  a  Gaussian  form  (as  a  function  of  wave  number,  k)  for  the 
angular  power  spectrum  as  well.  I f  we  now  characterize  the  spatial 


» 
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decorrelation  function  in  a  particular  direction  by  the  value  of  A, 

denoted  A  ,  for  which  it  has  decreased  to  1/e  [i.e.,  D(A  )  =  2],  the 
e  e 

half  width  of  the  angular  spectrum  in  this  direction,  0  ,  again 

e 

characterized  by  the  1/e  point,  is  given  by 

9  =  arc  sin  — — —  ,  (36) 

S  nf/A~ 

e 

where  c  is  the  speed  of  light  and  f  is  the  frequency.  The  apparent 
complexity  of  Eq.  (36)  is  a  consequence  of  the  width  being  expressed  in 
terms  of  the  angular  variable,  0,  rather  than  a  component  of  the 
wave  number,  k. 


SECTION  6 


TEMPORAL  DECORRELATION 

Temporal  decorrelation  of  the  signal  is  caused  by  the  evolution  of 
the  irregularities  through  which  the  signal  propagates.  It  can  be 
evaluated  in  a  manner  similar  to  that  used  to  determine  the  spatial 
decorrelation,  discussed  in  the  previous  section.  For  temporal 
decorrelation,  however,  the  differences  between  the  media  encountered 
along  the  same  path  at  two  times,  rather  than  those  found  between  two 
paths  at  the  same  time,  are  of  interest. 

If  a  fixed  spatial-irregularity  structure  is,  as  we  assume, 
essentially  convected  across  the  path  during  the  time  necessary  for  the 
signal  to  decorrelate,  the  effect  is  the  same  as  if  the  path  were  being 
shifted  spatially.  In  this  case,  the  calculation  of  the  temporal 
decorrelation  becomes  very  similar  to  that  described  previously  for 
spatial  decorrelation.  To  apply  the  calculation  to  temporal 
decorrelation,  it  is  necessary  only  to  modify  the  details  about  how  the 
path  is  to  be  displaced  in  a  suitable  manner. 

The  simplest  case  is  that  of  a  uniform  convection  velocity  at  all 
points  along  the  path.  In  this  case,  a  uniform  displacement  of  the 
entire  path  in  the  direction  of  the  velocity  vector,  with  the  magnitude 
of  the  displacement  proportional  to  time,  is  appropriate.  The 
modification  of  the  phase-structure-function  calculation  required  to 
produce  such  a  displacement  is  readily  achieved  by  redefining  the 
displacement  parameter,  y,  of  Eq.  (23).  This  parameter  is  considered 
proportional  to  a  time  separation.  At: 


The  proportionality  constant  v  constitutes  aR  effective  speed  for 

drift  of  the  irregularities  across  the  path.  For  a  drift  in  an 

arbitrary  direction,  v  =  (v.,v  v  ),  v  c,  is  given  by 

o  <p  r  erf 

2  2 
Cv  -Bv  v  +  Av 
~  s ,  s  ,  s  s 

V  ff  -  ♦  \ - *  .  (38) 

AC  -  B  /4 


which  corresponds  roughly  to  Eq.  (23).  In  Eq.  (38), 


v  =  v,  -  tan 
s  i 


6\  V9 


with  v  =  (v  ,v  ),  much  as  in  Eq.  (28).  The  vector  a  in  Eq.  (39)  is  a 
l  <p  r  t 

unit  vector  in  the  r—  4»  plane  that  is  directed  along  the  projection  of 

the  wave-normal  vector,  k,  into  this  plane,  8  being  the  angle  between  k 

and  the  perpendicular,  directed  along  i,,  to  this  plane.  Thus, 

-  9 

a,  =  ±i  ,  depending  upon  whether  the  portion  of  the  ray  going  up  (+) 
or  going  down  (-)  is  being  considered. 

The  temporal  variation  of  the  medium  can  also  be  discussed  in  terms 
of  the  distortion  of  the  signal  frequency  spectrum  (i.e.,  Doppler 
spread)  that  it  produces.  As  in  the  case  of  the  spatial-correlation 
function  and  angular  power  spectrum  discussed  in  Section  5,  the  Fourier- 
transform  relationship  that  connects  the  signal  power  spectrum  and  the 
time  correlation  function  yields  particularly  simple  results  when  the 
power-law  index,  v,  of  the  ionospheric  irregularities  [Eq.  (2)]  has  a 
value  of  1.5.  Again,  the  form  of  both  functions  becomes  Gaussian  in 
this  circumstance,  and  the  1/e  decorrelation  time,  t,  and  the  spectral 
1/e  half  width,  f^,  are  related  in  essentially  the  same  manner  as  are 
their  spatial  counterparts: 


SECTION  7 


FREQUENCY  DECORRELATION 

Decorrelation  of  the  signal  frequency  spectrum  and  a  corresponding 
time  (delay)  spread  of  the  signal  waveform  also  result  from  scatter  by 
ionospheric  irregularities.  Several  effects  can  be  identified  that  may 
contribute  to  the  delay  spread.  Two  of  these  are  familiar  from  the 
theory  previously  developed  for  unrefracted  line-of-sight  paths.  The 
third,  which  bears  a  close  relationship  to  the  frequency  spread  and 
angular  spread  that  have  already  been  discussed  in  previous  sections  of 
this  report,  is  potentially  important  only  for  refracted  ray  paths. 
Consequently,  it  appears  not  to  have  been  considered  previously. 

We  first  note  that  time  spread  can  result  from  the  difference  in 
travel  time,  or  group  delay,  between  the  unscattered  ray  path  and  those 
that  join  the  transmitter  and  receiver  to  scattering  points  lying  to  one 
side  of  the  unscattered  path.  The  magnitude  of  this  component  of  the 
time  spread  is  determined  largely  by  the  difference  in  geometric  path 
length  between  the  two  paths. 

Considerable  simplification  in  calculating  this  delay-spread 
component  can  be  achieved  by  observing  that  the  constraints  imposed  by 
refraction  on  vertical  displacement  of  the  ray  path  severely  limit  the 
sensitivity  of  the  signal  to  scatter  by  the  vertical  component  of  any 
irregularities  present.  This  insensitivity  is  further  enhanced  by 
irregularity  structures,  such  as  those  encountered  in  the  ionosphere  at 
mid  and  high  latitudes,  that  have  much  larger  scales  vertically  than 
horizontally.  Thus,  to  good  approximation,  we  need  consider  only  the 
lateral  spread  in  the  signal  angular  spectrum  to  determine  this 
component  of  the  delay  spread.  If  we  again  assume  that  the  significant 
scatterers  are  concentrated  near  midpath,  a  simple  geometric  estimate, 
ignoring  earth  curvature,  of  the  delay  spread,  tg,  yields 


(41) 


t 

s 


D  2- 
t—  tan  0 
2c  e 


where  D  is  the  path  length,  and  0^  is  the  e  *  point  on  the  lateral 
spread  owing  to  scatter  of  the  signal  angular  spectrum,  which  was 
discussed  in  Section  5. 


We  have  assumed  (for  the  present)  that  the  major  contribution  to 
the  delay  spread  will  generally  result  from  the  effect  just  discussed. 
This  expectation,  however,  remains  to  be  placed  on  a  firm  footing. 

Thus,  it  seems  appropriate  to  discuss  briefly  the  other  possible  sources 
of  delay  spread  and  how  they  might  be  assessed. 

Delay  spread  can  also  result  from  random  fluctuations  of  the 
integrated  electron  density  along  the  unscattered  ray  path.  Such 
fluctuations  on  line-of-sight  paths  cause  a  corresponding  fluctuation  in 
the  group  delay  that  constitutes  a  delay  spread.  For  refracted  ray 
paths,  however,  the  effect  of  such  fluctuations  separate  less  clearly 
from  those  associated  with  scatter  from  points  to  one  side  of  the  path, 
discussed  above.  It  remains  to  be  determined  if  a  simple  integration  of 
the  mean-square  refractive-index  fluctuations  along  the  unscattered  ray 
path  will  adequately  describe  this  contribution  to  the  delay  spread. 

A  third  component  of  the  delay  spread  is  most  easily  addressed  in 
the  frequency  domain,  as  a  frequency  decorrelation.  In  contrast  to  the 
situation  at  higher  frequencies,  in  which  refraction  is  generally 
considered  to  be  negligible,  a  shift  in  signal  frequency  at  HF  displaces 
the  path  to  first  order.  This  spatial  separation  implies  a  frequency 
decorrelation  in  consequence  of  the  random  spatial  fluctuations  of  the 
medium  through  which  the  signal  propagates. 

This  decorrelation  can  be  evaluated  in  much  the  same  way  as  is  that 
from  a  spatial  displacement  of  the  receiver,  namely,  by  calculating  the 
phase-structure  function,  as  discussed  in  Sections  4  and  5.  The  spatial 
separation  produced  at  each  point  along  the  unperturbed  ray  path  by  a 
small  frequency  displacement  is  first  determined.  The  incremental 
contribution  to  the  phase-structure  function  from  this  displacement  is 


1  ') 


SECTION  8 


EFFECTS  IN  A  NOMINAL  ENVIRONMENT 


0 

Sample  calculations  made  in  the  course  of  implementing  the  model 
for  use  in  conjunction  with  the  HF  Channel  Probe  experiment  serve  to 
illustrate  the  model.  The  results  of  some  of  these  calculations  are 
presented  here.  The  values  assumed  for  the  model  parameters  in  these 
calculations  are  given  in  Table  2.  These  values  are  representative  of 

Table  2.  Parameter  values  for  a 
sample  environment. 


Parameter 

Value 

V 

1.5 

qo 

10"4  m_1 

a" 

10 

b' 

1 

<an2>/n2 

e  e 

0.8 

Veff 

400  mis 

either  the  natural  environment  or  one  disturbed  by  a  high-altitude 
nuclear  detonation.  The  value  of  0.8  assumed  for  the  fluctuation 
variance  is,  however,  probably  excessively  large  for  the  natural 
environment.  Measurements  using  the  DNA  HILAT  satellite  have,  for 
example,  shown  0.2  to  be  a  more  representative  value  for  this  parameter 


in  the  polar-cap  ionosphere,  which  is  a  relatively  severe  example  of 
naturally  produced  structure.  The  400  m/s  drift  speed,  on  the  other 
hand,  is  low  for  the  nuclear  environment;  1000  m/s  might  be  more 
appropriate  as  a  nominal  value  for  the  latter. 

A  number  of  calculations  were  made  for  a  propagation  path  1918  km 
in  length  (the  length  of  the  path  instrumented  experimentally)  using  the 
two  profile  models  shown  in  Figure  2.  Both  models  have  vertical- 


PROFILE  A 


PROFILE  B 


Figure  2.  Model  electron  density  profiles. 


C.  Kino,  personal  communication. 
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Incidence  critical  frequencies  of  4.74  MHz  at  a  peak  height  of  310  km. 
Profile  A  has  a  semithickness  of  91  km;  Profile  B,  19!  km.  The  path 
also  was  taken  to  lie  in  the  geomagnetic  meridian  plane,  as  was 
appropriate  for  the  intended  experiment,  with  a  magnetic-dip  angle  of 
80.43  deg.  The  spatial  and  temporal  correlations  and  the  related 
quantities  calculated  for  a  signal  propagated  over  this  path  Cat  a 
frequency  of  10.12  MHz)  are  shown  in  Table  3.  The  1/e  distances  given 


Table  3.  Calculated  scatter  characteristics 
for  a  sample  environment 


Profile 


Doppler 

1/e  Time  Spread 
(ms)  (Hz) 


1/e  Distance  (m)  |  1/e  Angle  (Deg) 


453.0  32.4  1.2 


624.0  58.9  0.9 


for  the  spatial  decorrelation  are  those  along  the  principal  axes  of  the 
spatial  correlation  ellipse;  for  a  path  in  the  meridian  plane,  these 
axes  lie  in  the  plane  containing  the  path  (II)  and  perpendicular  to  this 
plane  (  ).  The  delay  spreads  given  in  the  last  column  of  the  table  were 
calculated  from  the  angular  spreads  using  Eq.  41.  The  spatial  and 
temporal  correlation  functions  are,  according  to  Eq.  35,  Gaussian  in 
these  examples  by  virtue  of  the  value  of  1.5  assumed  for  the  power-law 
index  of  the  spatial  irregularity  spectrum,  and  the  linearizing 
assumptions  made  regarding  the  path  displacement  in  the  derivation  of 
Eq.  33. 

As  expected,  given  the  magnetic  field  orientation  and  the 
elongation  of  the  irregularities  along  this  field,  the  angular  spreads 
shown  in  Table  3  are  much  larger  horizontally  than  they  are  in  the  plane 
of  the  path.  It  is  also  interesting  that  Profile  B  gives  smaller 
spreads  in  both  angle  and  Doppler  than  does  Profile  A.  This  behavior  is 
consistent  with  the  ray  path  penetrating  to  somewhat  higher  electron 
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